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(54) Methods and apparatus for detecting certain compounds 

(57) The detection of compounds containing both nuclei of a kind which exhibit nuclear magnetic resonance (nmr) and 
nuclei of a kind which exhibit nuclear quadrupole resonance (nqr), for example narcotics and explosives containing both 1 H 
and "N nuclei, is carried out by subjecting a substance or object suspected of containing such a compound to a magnetic 
field and repetitively to three parallel activities A,H and B. Activity A comprises two similar rf pulse sequences P A and P B at 
the nmr frequency. Activity H is a temporary change in the value o1 the magnetic field, which may be a reduction to zero. 
Activity B is the application of additional rf signals. Activity H and/or activity B are different during each of the respective 
pulse sequences P A and P r The activity B rf signals are at one or more frequencies selected so that in combination with 
the activity H in at least one of the pulse sequences the subsequent nmr signal is modified H such a compound is present 
due to energy coupling between the two kinds of nuclei. The difference in magnitude of the nmr signals obtained following 
the two types of rf pulse sequence is then measured. 
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129581 

METHODS AND APPARATUS 
FOR DETECTING CERTAIN COMPOUNDS 
The invention relates to methods and apparatus for detecting 
certain compounds. More particularly the invention relates to 
methods and apparatus for detecting compounds containing both 
nuclei of a. first kind which exhibit nuclear magnetic resonance 
05 (nmr) and nuclei of a second kind which exhibit nuclear 
quadrupole resonance (nqr). Examples of such compounds are a 
range of narcotics and explosives such as cocaine, heroin, TNT, 
RDX and PETN. 

Nuclei of the first kind usually have a spin quantum number 

10 of 172 and include the nuclei. 1 H, 19 F, 31 P which give rise to 
strong and easily detectable nmr responses in the presence of a 
magnetic field. The nmr frequency of such nuclei is given by 
u=yH where co is the angular frequency, y the gyromagnetic 
ratio specific to that nucleus and H Is the intensity of the 

15 magnetic field. The nmr frequency thus depends on the value of 
the applied magnetic field. Nuclei of the second kind have a 
spin quantum number greater than 1/2 and include 11 B, 14 N, 
27 A1, 35 C1. They exhibit characteristic frequencies known as 
their quadrupole resonance frequencies. These frequencies are 

20 usually much less dependent on the intensity of the applied 
magnetic field. For a given nucleus these frequencies depend on 
the molecular structure within which it is contained and 
therefore provide -a sensitive identification of individual 
compounds containing such nuclei. The examples of narcotic and 

25 explosive substances quoted above are compounds containing both 
the X H nucleus and the 14 N nucleus and the invention is 
particularly useful in the detection of compounds containing 
these nuclei, although the Invention can also be used to detect 
compounds having different combinations of nuclei. 

30 It has been suggested that by causing any nqr response, if 

present, to change the Intensity of an induced nmr signal by 
transfer of energy between the two kinds of nuclei the presence 



of substances containing both kinds of nuclei could be detected. 
This phenomenon is known as double resonance and can be made to 
occur, for example, when the nmr frequency is made to coincide 
with the nqr frequency, a condition known as level-crossing. 

05 One such suggestion is described in British Patent 

No. 2057135. In that specification it is proposed to generate an 
nmr echo signal by the application of a pair of rf pulses of 
appropriate frequency in a magnetic field. The magnitude of this 
echo signal is measured. The magnitude of the magnetic field is 

10 then changed to a value at which the nmr frequency coincides with 
the nqr frequency of a nucleus of the second kind in a compound 
which it is desired to detect so shortening the nmr spin-lattice 
relaxation time, a phenomenon known as cross-relaxation. The 
magnetic field intensity is then returned to its initial value 

15 and a similar pair of pulses applied and the echo signal again 
measured. If there has been cross relaxation when the intensity 
of the magnetic field was at its changed value then the echo 
signal response after the second pair of pulses will be 
different. This will indicate the presence of the compound. 

20 However, the two-pulse sequence described is intended to generate 
a so called "solid-echo" signal and to be effective the two 
pulses must be in phase quadrature. 

Another method is described in US Patent No. 4,296,378. 
However the method described there does not distinguish between 

25 substances in which cross-relaxation occurs and those substances 
containing nuclei of the first kind with similar nmr spin-lattice 
relaxation tiroes. 

Another suggestion for detecting the presence of certain 
substances is described in published PCT Application 

30 No. W084/04173. The method described there is based on the 
assumption that in substances of interest the *H nuclei have long 
spin-lattice relaxation times (l x ) and short spin-spin relaxation 
times (T 2 ). There are, however, many hazardous materials which 
it may be desired to detect which do not fit into this category. 

35 One example is ammonium nitrate in which T x is relatively short 



and T 2 is relatively long. Another example is cocaine where T t 
is short both in high and low magnetic fields. 

It is an object of the invention to provide a method and 
apparatus capable of detecting any one of a large range of 

05 compounds in which double resonance can be made to occur. 

According to the invention in one aspect a method of 
detecting the presence of one or more compounds containing both 
nuclei of a first kind which exhibit nmr and nuclei of a second 
kind which exhibit nqr comprises subjecting material suspected of 

10 containing such a compound to a magnetic field and successively 
and repetitively to two types of sequence, each sequence 
comprising an rf pulse or pulse train at the nmr frequency of the 
first kind of nucleus in the said magnetic field followed by a 
time interval in which in at least one of the sequences there is 

1.5 a predetermined change in the magnitude of the said magnetic 
field, and during the time interval in at least one of the types 
of sequence applying rf signals at one or more frequencies 
selected so that in combination with the aforesaid change in the 
magnitude of the magnetic field the subsequent nmr signal is 

20 modified if such a compound is present, and measuring the 
difference in the magnitude of the nmr signals obtained following 
the two types of sequence. 

Where rf pulse trains are employed, then these trains can 
comprise spin, locking pulse sequences comprising 90° pulses of 

25 the form 90* o -T«(90 p 9O -2T) n where t is a time delay and n is an 
integer of value one or more. Alternatively the rf pulses may 
comprise single or multiple 90° pulses. The two types of 
sequence may be applied alternately. Alternatively, a plurality 
of one type of sequence is first applied and is then followed by 

30 an equal plurality of the other type of sequence. 

In one embodiment of the invention the changes in the 
magnitude of the magnetic field are changes to a zero or low 
value and rf signals at the nqr frequency of one or more 
compounds which it is desired to detect are applied in the time 

35 interval of one of the types of sequence. In addition rf signals 



not at the nqr frequency of any of the compounds of interest may 
be applied during the corresponding time interval of the other of 
the types of sequence. 

In another embodiment of the invention the change in the 

05 magnitude of the static magnetic field during one of -the types of 
sequence is a stepped change with at least one of the steps being 
held at a value in which level crossing between nuclei of the two 
kinds takes place if corresponding ones of said compounds are 
present. The rf signals are applied when the magnetic field 

10 reaches a zero value in order to destroy any residual nmr signal. 

In yet another embodiment of the invention the changes in the 
magnitude of the static magnetic field are different in the 
respective time intervals of each of the types of sequence. One 
of the changes is such as to cause the nmr frequency of the first 

15 kind of nuclei to coincide with one of the nqr frequencies of 
nuclei of the second kind in a compound which it is desired to 
detect. Simultaneously therewith rf signals are applied at a 
frequency corresponding to a different nqr frequency of such 
nuclei . 

20 According to the invention in another aspect apparatus for 

detecting the presence of one or more compounds containing both 
nuclei of a first kind which exhibit nmr and nuclei of a second 
kind which exhibit nqr comprises: 

means for applying a magnetic field to material to be tested; 
25 means for applying rf pulses to such material of predeter- 

mined radio frequency, intensity, duration and phase at the nmr 
frequency of nuclei of the first kind; 

means for temporarily changing the magnitude of the said 
magnetic field in a predetermined manner; 
30 means for applying rf signals to such materials during the 

time of said temporary change of magnitude of the said magnetic 
field said signals being of one or more predetermined 
frequencies; and 

means for measuring the free induction signal from nuclei of 
35 the first kind. 



In carrying out the invention separate rf coil means may be 
provided to which the said rf pulses and the said rf signals 
respectively are applied. Two separate receiving channels may 
also be provided to receive different ones of free induction 
05 signals from nuclei of the first kind. 

In order that the invention may be more fully understood 
reference will now be made to the accompanying drawings in which:- 
Figure 1 shows a timing diagram of rf pulses, switched 
magnetic fields and applied rf signals that can be used in a 
10 method embodying the invention, 

Figure 2 shows in detail one rf pulse train that may be used 
in the method of Figure 1, 

Figures 3,4 and 5 show alternative timing diagrams of pulses, 
switched magnetic, fields and applied rf signals used in methods 
15 embodying the invention. 

Figure 6 shows apparatus embodying the invention in block 
diagrammatic form, and 

Figures 7, 8a, 8b, 9 and 10 give results obtained by using the 
apparatus of Figure 5 in accordance with one or other of the 
20 techniques shown in Figures 1,3,4 and 5. 

•Referring now to Figure 1, there are three parallel 
activities that take place, which are- labelled A, H and B. The A 
activity comprises the application of two types of sequence 
labelled Pa and Pg. Both types of sequence include similar rf 
25 pulse trains indicated by the composite pulse envelope P$l 
followed by a time interval x'. The activity H relates to the 
application of a magnetic field which is normally held at a 
constant value H 0 but within the interval x' it is switched to a 
zero value for a time x' p which is less than x'. There is a 
30 delay time x<j between the end of interval x' p and of the interval 
x 1 and a further delay time given by x'-x^-x'p between the end of 
pulse train P$l and the commencement of interval x'p. In the 
time Interval x* p an activity labelled B takes place. Activity B 
consists of the application of an rf signal as shown by the B 
35 line. This signal may have one or more different frequencies as 
will be described hereafter. 
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The pulse envelope P$l itself comprises single or multiple 
90° pulses or one or more- pulse trains in the form of a multiple 
pulse spin-locking sequence which is shown in more detail in 
Figure 2. The spin-locking pulse sequence comprises a first 90° 

05 rf pulse followed by a train of further 90° pulses whose rf phase 
is in quadrature with the first pulse. The first pulse is 
written as 90° o and subsequent pulses as 90% o . The time spacing 
between the first pulse and the commencement of the series of 
90° 9O pulses is x and the time spacing between each of the 90° 9O 

10 pulses of the series is 2x. The complete pulse sequence can 
therefore be written as 90 0 0 -T-(90° 90 -2T> n where n is an integer 
of value one or more. Provided that the frequency of the rf 
pulses is chosen relative to the magnitude of the applied 
magnetic field H 0 so as to satisfy the relationship w=y h o where w 

15 is the angular frequency of the rf pulses and y is the 
gyromagnetic ratio of nuclei of the first kind, then nuclear 
magnetic resonance will- be induced and a free induction signal 
will be generated after each individual pulse of the sequence. 
This signal peaks during the time interval 2x between successive 

20 90° pulses and its peak value is measured and stored for later 
use. The entire pulse train that Is shown in Figure 2 is 
designated by the pulse envelope P$l 1° Figure 1. 

Turning now back to Figure 1 attention is directed to the 
line marked B which shows the application of rf signals during 

25 the time that the magnetic field H has a zero value. In the time 
interval in the sequence P^ this signal is chosen to have a 
frequency which is not equal to any of the nqr frequencies of 
compounds suspected of being included in the material under 
investigation. However, in the sequence Pg one or more 

30 frequencies are applied either simultaneously or sequentially 
which are made equal to at least one of the nqr frequencies of 
suspected compounds. 

In use of the timing diagram shown in Figure 1 the sequences 
Pa and Pb are applied successively and repetitively. The actual 

35 overall sequence can for example comprise a series of Pa 



sequences followed by an equal series of Pg sequences, as 

p A p A p A- • - P B P B P B Alternatively the sequences can be applied 

alternately i.e. p A p B p A p B--- 

Measurements of the free induction signals for both the Pa 

05 and the Pb sequences are taken using quadrature phase-sensitive 
detectors and are added separately for each sequence. The summed 
free induction signals may be Fourier transformed to obtain the 
absorption-mode signal which is then integrated to improve the 
signal-to-noise ratio. If no suspected material is present then 

10 there will be no difference in the magnitudes of the two 
integrations. However, if any such material is present then 
there will be a difference. The existence of this difference 
indicates the presence of one or more suspected compounds 
containing nuclei at the frequencies of the B signals applied 

15 during pulse sequence Pg. 

To facilitate understanding of the way in which the above 
timing diagram functions, the two kinds of nucleus that may be 
present in a material being examined are termed in the following 
discussion as A nuclei and B nuclei. The application of 

20 successive and repetitive sequences P/\ and Pg may generate 
pseudo-equilibrium conditions in the A-spin system in the sense 
that the A populations are maintained away from their thermal 
equilibrium values. Use of pseudo-equilibrium conditions has the 
advantage that the material being examined can be immediately 

25 subject to the above described pulse sequences and there is no 
need to wait for the material to acquire a finite magnetisation 
tn the applied magnetic field. 

The magnetic field as shown in the H line in Figure 1 is 
cycled to zero adiabatical ly but slowly enough to ensure that the 

30 population ratios of both the A and B spin systems become equal 
at level crossing during the change. Level crossing will occur 
as the magnetically-split A-spin levels, equal to y^H* collapse 
to zero (or increase from zero) and cross over the B-spin levels, 
the quadrupole interaction of which ensures a non-zero frequency, 

35 cog , even in zero field. At level crossing, when the condition 



- 8 - 

"B =YA H * s satisfied, the A and B spin temperatures Ta,Tb defined 
by the population ratios 



N A upper 



=e - E A /kT A. 



N B upper 



r E B'" Bt 



N A lower 



N B lower 



become equal in a time T^p, the cross-polarisation time. The 
field cycling must therefore be slow enough in this train to. 
05 ensure that the time spent at the level-crossing condition be 
longer than Trjp. 

During the time interval x p ' in sequence Pa, off-resonant rf 
radiation of angular frequency w, with or without frequency 
modulation may be applied, as shown by the B line in Figure 1, to 

10 the B spins in zero or low magnetic field. Frequency w is 
off-resonant with respect to the B quadrupole resonance 
frequencies wq. During time x p ' in sequence Pg the B-spins are 
subject, in zero or low magnetic field, to resonant rf radiation 
with or without frequency modulation at saturation level at one 

15 or more of their quadrupole resonance frequencies uq (1 > ,uq (2) , at 
the temperature of the object. For a three-level spin system 
like 14 N (1=1), with three resonant frequencies labelled 
o x ,wy f w z , usually w x and/or u z will be saturated and the power 
spectrum of the B pulse can contain any number of different oj x 

20 (and/or u 2 > frequencies characteristic of all nitrogen-containing 
compounds in the sample whose presence is to be detected. In 
this way all such compounds are simultaneously irradiated. On 
remagnetising adiabatically the saturated i4 N (B) levels pass on 
their energy to the l H (A) levels on level crossing. The 

25 magnetic field and irradiation cycle may be repeated several 
times during the time t p * to produce multiple level-crossing and 
enhance sensitivity. Finally, the reduced X H signal is observed 
by a free induction signal accumulator in the time sequence Pg 
relative to the corresponding signal in the time sequence Pa- 
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This timing sequence shown in Figure 1 is related to that 
used in double resonance by level crossing experiments (drlc). 
The technique is. particularly suitable for materials in which 
both the A nuclei and the B nuclei have medium to long 

05 spin-lattice relaxation times T x and with a reasonably high 
proportion of B nuclei to A nuclei. This occurs for the A 
nucleus l H and the B nucleus 14 N in many explosives. 

For compounds in which the cross-polarisation time T^p is 
long but the B nuclei are relatively abundant , as in the 14 N0 2 

10 group in some explosive^, the technique of Figure 3 may be used. 
In this technique the pulse train P$|_ is followed by time 
intervals in which the magnetic field is changed as in Figure 1, 
but the nature of the change and the nature of the rf signals 

15 applied in the time intervals are different. For the highest 
sensitivity time x r should be equal to or longer than the time l x 
for the A nucleus. The magnetic field cycle is adiabatic but may 
be too fast to produce equalisation of the A and B spin 
temperatures during level crossing. During time x p ' in sequence 

20 Pa the magnetic field is cycled straight to zero and 
simultaneously the A dipolar system saturated with low frequency 
radiation, typically 10 kHz, to destroy any remaining A 
magnetisation. During time x p \ in sequence Pg the same sequence 
of events is used except that the field cycle is now conducted 

25 slowly enough for equalisation of the A and B spin temperatures 
at level crossing to occur e.g. by halting for a time Tcp or 
longer at one or more level crossing fields in which the A nmr 
frequency matches in turn one or more of the B quadrupole 
resonance frequencies of each of the materials to be detected, 

30 both during demagnetisation and remagnetisation. During the 
latter the B nuclei quadrupole levels share their surplus 
populations, acquired during the former, with the A nuclei spins 
in sequence, producing an increase in the A signal. If T x (A) is 
short, time x" (and xj) may be made so short as to reduce or 

35 eliminate many A responses from solids usually encountered in the 
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objects to be examined; only materials containing B nuclei whose 
quadrupole resonance frequencies are matched at the level 
crossing steps will then give significant signals. 

For compounds in which T 1 (A)>>T 1 (B) t x- can be made much 

05 shorter than T 1 (A) and the magnetic field steps, now needed on 
remagnetisation only, are made to last much longer than the 
cross-polarisation time T^p so that considerable re-polarisation 
of the A spins by relaxation recovery occurs in the 
level-crossing field due to thermal contact with the B spins. 

10 Away from the level crossing condition re-polar1sation is much 
slower. Again, fast repetition of the pulse units may reduce the 
background signal considerably. 

These pulse trains may be called double resonance by 
adiabatic remagnetisation (dram). 

15 The technique of Figure 4 is appropriate for compounds in 

which the time T, for the A nucleus 1s much longer than T x for 
the B nucleus and the B nucleus is in low abundance, as for 14 N 
in heroin and other narcotics. Under these circumstances, if the 
A Zeeman levels are matched in low magnetic fields with a pair of 

20 quadrupole-split B levels, after a short period T^p 1n which the 
two spin systems come to a common spin temperature, the combined 
system A+B thereafter relaxes with a spin-lattice relaxation time 
(T 1 ) e ff which is the weighted mean of the two individual 
relaxation times T X (A) and T 1 (B) 



<T 1 ) ef f= 



<-£->JL + (-i-)O- 



-1 



Uc T](B) 1+c T 1 <A>J 

25 in which c is equal to (or proportional to) the number ratio of B 
to A spins Nb/Na- If c T 1 (A)/T 1 (B)>>1 , the method has high 
sensitivity for the detection of B spins. 

In this sequence as shown in Figure 4, t' is set close to J x 
in the field H. During time tp' in sequence Pa the magnetic 

30 field is cycled adiabatically to a magnetic field in which the A 



magnetic resonance frequency is close to, but not equal to, one 
of the B quadrupole resonance frequencies characteristic of the 
material it is wished to detect. During time xp* in sequence Pb 
the magnetic field is cycled to a value such that the A magnetic 

05 resonance frequency, is exactly matched to the peak- B frequency 
response. The A-nucleus spin-lattice relaxation is much faster 
during sequence Pb than during sequence Pa, resulting in a faster 
loss of polarisation during Pb than Pa and very different 
accumulated signals at the two detectors. This difference may be 

10 even further increased by level crossing effects in the 
appropriate circumstances e.g. when T A (A) is short. 

In a simpler version of this cycle which is applicable when 
the measurement field H 0 and the cross-relaxation field H are 
close in magnitude, pulse envelope P$l comprises a saturation 

15 recovery sequence of- 90° pulses with separations comparable to 
T X (A) in field H 0 . In sequence Pa, the magnetic field is set to 
H 0 and during Pb to H for an equal time, the A nmr frequency 
being changed accordingly. The A-nucleus spin-lattice relaxation 
is much faster during sequence Pb than during sequence Pa, 

20 resulting in a partial recovery of the A signal during sequence 
Pb relative to that in Pa and different accumulated signals at 
the two detectors. 

In another version of the cycle, the magnetic field during 
time Tp' of sequence Pb is set to half the magnetic field 

25 required for matching so as to observe the narrower two-proton 
relaxation jumps. 

In all versions of this cycle appropriate to spin-1 B nuclei 
such as 14 N, rf signals are applied in time Tp 1 in both sequences 
Pa and Pb of frequency equal to another quadrupole resonance 

30 frequency differing from that which is being cross-relaxed. If, 
for example, the cross-relaxation magnetic field is matched to 
one nqr frequency «b' irradiation at another nqr frequency «b" 
. causes positive enhancement of the recovered A magnetisation, 
providing thereby both an increase in sensitivity and a highly 

35 sensitive identification of the material present in the object. 



This technique Is related to a group of methods known 
collectively as double resonance by cross relaxation (drcr); it 
can detect 14 N signals in crystalline solids, in plastics and 
also in liquids and. gels with long correlation times OCT 5 to 
05 10- 6 s>. 

The technique of Figure 5 is appropriate to compounds in 
which the quadrupolar nucleus B has half-integral spin (eg I« 3 /2> 
and T^B) is short, such as the 35 C1, 37 C1 nuclei in cocaine 
hydrochloride. During time x p ' in sequence Pa as shown in 

10 Figure 5 the magnetic field is switched as fast as possible down 
to or close to zero and intense rf radiation applied at the B 
quadrupole resonance frequency (wq) of the compound to be 
detected at the given temperature. During time x p ' in sequence 
Pb the cycle is the same except that the rf radiation is 

15 off-resonance by a carefully specified shift a=w-o>q. During the 
second cycle the effective frequency "seen" by the B nuclei in 
the rotating frame is.«-uQ=A, and if this lies within the A 
dipolar absorption edge the A spins become polarised provided 
that T^B) is short. Fast return to higher field with the 

20 minimum value of t<j "captures" • this polarisation within the 
A Zeeman levels resulting in an enhanced A signal if u«oq and a 
diminished (or inverted) A signal if w>uq. If the timing period 
x' is made so short that the A signal is saturated, any remaining 
signal in zero field being destroyed by an initial burst of 

25 radiation within the A dipolar absorption edge, the A detector 
will register no response as will the B detector unless the 
B-nucleus whose frequency is irradiated is present. As in the 
technique of Figure 3 it Is possible in this way to eliminate or 
reduce much of the background A-signal from other A-containing 

30 materials in the object under examination. 

This technique is related to experiments based on double 
resonance by thermal mixing (drtm). 

Apparatus suitable for carrying out the methods shown in 
Figures 1,3,4 and 5 is shown in Figure 6. The apparatus has two 

35 rf channels. These channels provide the pulses and the signals 



for the A activity and the B activity respectively. The A rf 
channel comprises an A frequency source 1a. a pulse programmer 
and gate 2a, an rf power amplifier 3a and a coupling network 4a 
which couples the output of the rf amplifier 3a to an rf coil 

05 5a- Likewise the B rf channel comprises a B frequency source 1b. 
a pulse programmer and gate 2b, an rf power amplifier 3b and a 
coupling network 4b which couples the output of power amplifier 
3b to an rf coil 5b- 

The two coils 5a and 5b are positioned to provide rf fields 

10 to at least part of material under investigation and in addition 
a dc magnet 6 of controllable intensity provides a static 
magnetic field in which the material is immersed. The intensity 
of the magnetic field is controlled by a magnetic field cycling 
control 7. Received rf signals are detected in coil 5a and the 

15 detected signals are fed to a two-channel receiver 8 and thence 
to either one of two transient signal recorders 9a and 9b- 
Control of the frequency . sources 1a and 1b» of the pulse 
programmers and gates 2a and 2b, the magnetic field cycling 
control 7 and the transient recorders 9a and 9b is carried out by 

20 means of a computer 10 which takes the signals recorded in the 
transient recorders 9a and 9b processes them and provides an 
output 11. 

The apparatus of Figure 6 is programmed by computer 10 to 
operate in accordance with any one of the methods shown in 

25 Figures 1,3,4 or 5. The transient recorders 9a and 9b are 

triggered to detect the whole free induction decay signal, or in 
the case of a spin-locking sequence, only at the echo maxima. In 
a solid these signals decay in a time which is closer to T x p than 
T 2 , giving a much greater time for signal sampling. The sampling 

30 of the received signal can be performed for example by pulse 
gated integration or by fast digitisation and accumulation in the 
memory of the computer followed by Fourier transformation and 
Integration of the absorption signal. It Is desirable to use 
such accumulation techniques since the materials to be detected 

35 may be considerably diluted as, for example, being present in 



luggage and their signals will therefore be reduced by the 
presence of inert or harmless hydrogen-containing materials whose * 
response is not changed when the selected perturbation is applied 
to the B spin system. In the case of a large object being 

05 examined, such as a piece of luggage, it may be -necessary to 
examine it section by section by means of a sequentially switched 
array of A and B rf coils in order to keep the power requirements 
for the two amplifiers 3/\ and 3g within reasonable limits. 

Examples of results obtained from the use of the apparatus of 

10 Figure 6 are described below with reference to Figures 7,8,9 
and 10. In each case the most appropriate one of the methods of 
Figures 1,3,4 or 5 is used. The A nuclei in each case is the X H 
nucleus. The B nucleus is 14 N as present in various narcotics 
and explosives or else the 35 C1 or 37 C1 nucleus present in 

15 cocaine hydrochloride. 

The technique of Figure 1 detects a change in the X H 
magnetisation due to heroin under the following conditions. In 
the A channel a single 90° pulse is used of pulse width 3 \i% and 
the X H free induction decay at 40.18 MHz from a 3g sample is 

20 integrated in a box-car integrator of gate width 2 \is and time 
constant 0.1s. The magnetic field cycle p A p A- • - P B P B- • - is set to 
x , =4s,xp , =0.6s, T(j=0.1s and the B channel generates an rf field 
of amplitude about 0.2 mT. A change in signal of about 35 mV is 
produced when the frequency in the B channel is set to 4.024 or 

25 4.063 MHz at 295K. Figure 7 shows the complete high-frequency 
14 N spectrum produced by stepping the frequency after each 
sequence. The line is a doublet of splitting 39 kHz and has an 
overall width 59 kHz. Under similar conditions, cocaine free 
base ("crack") gives signals at 4.245 and 4.260 MHz of width 

30 30 kHz at the same temperature. The invention is therefore 
capable of detecting and distinguishing between these two 
compounds at room temperature. 

The techniques of Figures 3 and 4 give strong 14 NH 4 + and 
14 N0 3 * signals from 3g samples of ammonium nitrate under the 

35 following conditions. The A-channel setting is as in the * 



previous paragraph. The magnetic field cycle is that described 
1n Figure 4 with a single 1 H 90° pulse and t*=5.5s, t p '=0.3s and 
T(j=0.7s, the field being stepped after each sequence. As shown 
in the spectrum of Figure 8a the NH 4 + signals are detected at 

05 level-crossing magnetic fields equivalent to l H frequencies of 
48.9, 64.8 and 116.8 kHz. With the parameters changed to give 
t'-8.6s, t p '=1.2s and as before then the equivalent X H 
frequencies to detect the 14 N0 3 ~ Ion are at 103.9, 449.1 and 
549.0 kHz. This is shown in the spectrum of Figure 8b. A 

10 two-proton relaxation jump for the latter appears near 212 kHz 
and 1s sharp and strong enough to be used in detection. Mixtures 
of (NH 4 ) 2 S0 4 and (NH 4 )N0 3 show signals from both, whose relative 
intensity can be used to analyse for their relative proportions. 
The techniques of Figures 1 and 5 will detect signals from 

15 35 C1, 37 C1 1n cocaine hydrochloride. The A channel settings are 
as before. The magnetic field cycle has t**5s, Tp r =0.6s, T(j«0.1s 
and the B channel has rf radiation whose frequency Is stepped 
from 1.7 to 2.8 MHz with an H/ot f1eld of about o. 15 mT. The 
35 C1 signals are seen at 2.547 MHz and the signals from 37 C1 at 

20 1.998 MHz at room temperature. The complete spectrum is shown in 
Figure 9. This 35 C1 frequency is unique to cocaine hydrochloride 
and is most unlikely to be found in any other materials 
encountered in luggage, parcels or the like. 

Figure 10 shows the Vy lines near 3,400 kHz in the 14 N 

25 spectrum of the high explosive RDX. The recording was made by 
the same method as was used for NH 4 N0 3 but with timing settings 
of x*=50s, tp'sls and -C(j=0.7s. 
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CLAIMS 

1. A method of detecting, the presence of one or more compounds « 
containing both nuclei of a first kind which exhibit nmr and 

nuclei of a second kind which exhibit nqr comprising subjecting 
material suspected of containing such a compound to a magnetic 

05 field and successively and repetitively to two types of sequence, 
each sequence comprising an rf pulse or pulse train at the nmr 
frequency of the first kind of nucleus in the said magnetic field 
followed by a time interval 1n which in at least one of the 
sequences there is a predetermined change in the magnitude of the 

10 said magnetic field, and during the time interval in at least one 
of the types of sequence applying rf signals at one or more 
frequencies selected so that in combination with the aforesaid 
change in the magnitude of the magnetic field the subsequent nmr 
signal is modified if such a compound is present, and measuring 

15 the difference in the magnitude of the nmr signals obtained 
following the two types of sequence. 

2. The method as claimed 1n Claim 1 in which the rf pulse train 
comprises a sequence of 90° pulses. 

3. The method as claimed in Claim 2 in which the pulse train 
20 comprises a spin locking pulse sequence comprising 90° pulses of 

the form 90° o -T-(90° 9O -2T) n where x is a time delay and n is an 
integer of value one or more. 

4. The method as claimed in Claim 1 in which the rf pulse is a 
single 90° pulse. 

25 5. The method as claimed in any one of the preceding claims in 
which the two types of sequence are applied alternately. 
6. The method as claimed in any one of Claims 1 to 4 in which a 
plurality of one type of sequence is first applied and is then 
followed by an equal plurality of the other type of sequence. 

30 7. The method as claimed in any one of the preceding claims in 
which there are predetermined changes in the magnitude of the 
said magnetic field in both types of sequence. 
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8. The method as claimed In any one of the preceding claims in 
which the changes in the magnitude of the magnetic field are 
changes to a zero or low value and rf signals at the nqr 
frequency of one or more compounds which it is desired to detect 

05 are applied in the time interval of one of the types of sequence. 

9. The method as claimed in Claim 8 in which rf signals which 
are not at the nqr frequency of any of the compounds of interest 
are applied during, the corresponding time interval of the other 
of the types of sequence. 

10 10. The method as claimed in any one of Claims 1 to 7 in which 
the change in the magnitude of the static magnetic field during 
one of the types of sequence is a stepped change with at least 
one of the steps being held at a value in which level crossing 
between nuclei of the two kinds takes place if corresponding ones 

15 of said compounds are present. 

11. The method as claimed in Claim 10 in which the rf signals are 
applied when the magnetic field reaches zero value in order to 
destroy any residual nmr signal. 

12. The method as claimed in any one of Claims 1 to 7 in which 
20 the changes in the magnitude of the static magnetic field are 

different in the respective time intervals of each of the types 
of sequence. 

13 The method as claimed in Claim 12 in which one of the changes 
is such as to cause the nmr frequency of the first kind of nuclei 

25 to coincide with one of the nqr frequencies of nuclei of the 
second kind in a compound which it is desired to detect and 
simultaneously rf signals are applied at a frequency 
corresponding to a different nqr. frequency of such nuclei. 
14. Apparatus for detecting the presence of one or more compounds 

30 containing both nuclei of a first kind which exhibit nmr and 
nuclei of a second kind which exhibit nqr comprising: 

means for applying a magnetic field to material to be tested; 
means for applying rf pulses to such material of predeter- 
mined radio-frequency, intensity, duration and phase at the nmr 

35 frequency of nuclei of the first kind; 
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means for temporarily changing the magnitude of the said 
magnetic field in a predetermined manner; 

means for applying rf signals to such materials during the 
time of said temporary change of magnitude of the said magnetic 
05 field said signals being of one or more predetermined 
frequencies; and 

means for measuring the free induction signal from nuclei of 
the first kind. 

15. Apparatus as claimed in Claim 14 in which separate rf coll 
10 means are provided to which the said rf pulses and the said rf 

signals respectively are applied. 

16. Apparatus as claimed in Claim 14 or Claim 15 in which two 
separate receiving channels are provided to receive different 
ones of free induction signals from nuclei of the first kind. 

15 17. A method of detecting the presence of one or more compounds 
containing both nuclei of a kind which exhibit nmr and nuclei of 
a kind which exhibit nqr substantially as described herein with 
reference to Figure 1 or Figure 3 or Figure 4 or Figure 5 of the 
accompanying drawings. 

20 18. Apparatus for detecting the presence of one or more compounds 
containing both nuclei of a kind which exhibit nmr and nuclei of 
a kind which exhibit nqr substantially as described herein with 
reference to Figure 6 of the accompanying drawing. 
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75. A computer program product comprising: 

a computer usable medium having computer readable code 
embodied therein for establishing an identifier for at least a portion of a 
digital sequence comprising: 

computer readable program code devices configured to cause a 
computer to effect performing a function on said at least a portion of 
said digital sequence to produce a probabilistically unique symbol 
therefore; 

computer readable program code devices configured to cause a 
computer to effect establishing a correspondence between said at least 
a portion of said digital sequence and said probabilistically unique 
symbol; and 

computer readable program code devices configured to cause a 
computer to effect utilizing said probabilistically unique symbol as said 
identifier. 

76. The computer program product of claim 75 wherein said 
identifier and said corresponding at least a portion of said digital 
sequence are maintained in at least one data list. 

77. The computer program product of claim 76 wherein at least 
a portion of said identifier is utilizable as a pointer to a location of said 
corresponding at least a portion of said digital sequence within said at 
least one data list. 

78. The computer program product of claim 75 wherein said at 
least a portion of said digital sequence comprises at least a portion of a 
data file and said identifier is uniquely related to a content of said at 
least a portion of said data file. 

79. The computer program product of claim 75 wherein said at 
least a portion of said digital sequence comprises at least a portion of a 
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data stream and said identifier is uniquely related to a content of said 
at least a portion of said data stream. 

80. The computer program product of claim 75 wherein said at 
least a portion of said digital sequence comprises at least a portion of 
an executable file and said identifier is uniquely related to a content of 
said at least a portion of said executable file. 

81 . The computer program product of claim 75 wherein said 
computer readable program code devices configured to cause a 
computer to effect performing a function is carried out by: 

computer readable program code devices configured to cause a 
computer to effect hashing said at least a portion of said digital 
sequence to produce said probabilistically unique symbol. 

82. The computer program product of claim 81 wherein said 
computer readable program code devices configured to cause a 
computer to effect hashing is carried out by means of an industry 
standard digest algorithm. 

83. The computer program product of claim 82 wherein said 
computer readable program code devices configured to cause a 
computer to effect hashing is carried out by means of one of an MD4, 
MD5, SHA or SHA-1 algorithm. 

84. A computing environment substantially as herein before 
described with reference to the accompanying drawings as appropriate. 

A m e thod for manag i ng data comprising: 

producing a probab ili stically un i qu e id e ntif ie r for a d i g i ta l 

s e qu e nce; and 

comparing sa i d probab ili stically un i que i d e nt i fi e r to a list of oth e r 

i d e ntifi e rs corr e sponding to oth e r dig i ta l s e qu e ncos. 



38 

\\\DE - 53206/6 - H1104B1 vl 



MARKED -UP COPY 



85. A method for establishing an identifier for at least a portion 
of a digital sequence substantially as herein before described with 
reference to the accompanying drawings. 

Th e method of claim 84 furth e r compr i o i ng: 
5 adding said probab i l i st i cal l y un i que i dont i fi e r to said li st if sa i d 

probabi l i s tica ll y uniqu e id e nt i f ie r i s not pr e v i ous l y i n said l i st. 

86. A computer program product substantially as herein before 
described with reference to the accompanying drawings. 
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